Abstract
Two independent lines of evidence have been presented to the working groups and SC-CAMLR that claim to demonstrate that fishery-driven localised depletion of krill around pygoscelid penguin colonies has had a deleterious effect on their performance traits and demographic trends, that are equivalent to the impacts of climate variation. One study utilises 30 years of penguin foraging and reproductive performance measurements collected at two colonies in the South Shetland Islands while the other uses demographic rate changes derived from a comprehensive dataset of penguin population count data across Subarea 48.1 matched against acoustic measurements of krill biomass and krill catches at the gSSMU scale (Watters et al., 2020). The second uses estimated population trajectories across a wide range of penguin breeding colonies alongisde krill catches within 30km (Krüger et al., 2021). Both studies then explore the synergistic relationships to measurements of broad-scale climactic variation (El Nino Southern Oscillation;ENSO, and the Southern Annular Mode; SAM). Herein we provide a preliminary assessment of the efficacy of both approaches in drawing conclusions, that are now being used at the Commission level, as representing sound scientific advice. We demonstrate that several underlying assumptions in Watters et al. 2020 are contrary to the published scientific literature, and when the model syntax is re-written to reflect this, predicted penguin performance against long term expected means are substantially different to those presented to CCAMLR. STUFF ON KRUGER. While our preliminary assessment focuses on potential issues, future work will centre on considering competitive interactions both at appropriate time and space scales between the fishery as well as between a range of krill dependent predators beyond just pygoscelid penguins.
Concerns over the potential impact of localised depletion of krill through concentrated fishing effort on krill-dependent predators has been a topic of debate within SC-CAMLR and its Working Groups for many years (REF). Recently, two studies have been presented that suggest that local harvesting rates can impact predator performance to the same degree as poor environmental conditions (Watters et al., 2020) and when poor climactic conditions are coupled to locally high harvest rates the synergistic impacts on predators are evident (Krüger et al., 2021).
While both studies attempt to tackle the same overall problem, they do so using very different methodologies. Watters et al. (2020) exploit a considerable dataset; a substantial multi-species time series of penguin performance indices (including those collected under CEMP) collected over three decades at two sites (Cape Shireff on Livingstone Island and Copacobana on King George Island, South Shetland Islands; Figure 1) and over a decade of summer acoustic surveys that cover the at-sea distributions of Chinstrap, gentoo and Adélie penguins. Drawing in monthly krill catch statistics from the C1 Catch and Effort dataset and climactic data (Oceanic Niño Index; ONI), the authors use a hierarchical analysis of variance approach to estimate the variance in performance indices as a function of Local Krill Biomass (LKB), Local Harvesting Rates (LHR; the ratio of krill catch to LKB) and ONI. In contrast, Krüger et al. (2021) utilise a broader range of penguin colonies across the same three species throughout the Antarctic Peninsula area, in combination with their respective abundance survey estimates from an open-source database (www.penguinmap.org). The authors calculate population trends for appropriate sites, and using the CCAMLR C1 Catch and Effort data to extract annual catch values within a 30km radius of each colony. Finally, Krüger et al. (2021) use the trajectory of the trend (positive; increase versus negative; decrease) as a response in a binomial generalised linear mixed model using the summed annual catch, a lagged trend of the Southern Annular Mode and the SOMETHING ON ONI ???? to determine the relative contributions of each predictor and their interactive effects on population abundance trends. Both studies draw similar conclusions; that local harvesting levels of krill impact predators, and the degree of impact can either be similar to that of poor environmental conditions or have a synergistic impact when high local harvesting coincides with poor conditions.
These conclusions have been propogated into Commission documentation supporting the reformulation of the D1MPA proposal (CCAMLR-39/BG/02) as well as into Commission discussions (CCAMLR-39, Para 5.48 & Para 5.51). However, while the two studies have moved from Working Papers of EMM into the realms of the peer-reviewed literature, we have some areas of concern about the structuring of these studies that we feel warrant raising. Our preliminary review raises concerns unique to each study but also common across both, and we structure our paper accordingly. Firstly, we review Watters et al. (2020) and Krüger et al. (2021) through the lens of some of the ecological assumptions made versus the available evidence pertaining to them. Within the constraints of the data and analytical methods that are available from the studies, we also quantify how rationalising these assumptions to the evidence available impacts the conclusions drawn. We then highlight some overarching concerns applicable to both papers.
A key motivation for the paper appears to be to highlight the mismatch between the areal scales of fisheries management and ecological the interactions between fishing extractions and dependent predators. To do this, the authors create two strata aligned with groups of SSMU (gSSMU); gSSMU #1 including those SSMU inside the Bransfield Strait (APBSE and APBSW and gSSMU #2 incorporating SSMU north of the South Shetlands, including Elephant Island (APDPE, APDPW and APEI) represented in Figure 1. These gSSMU cover \(15,500nm^2\) and \(20,600nm^2\), respectively, and are used to characterise both krill biomass and harvesting rates that are “local” to the penguin colonies for which performance data are used. The reasoning behind scaling to gSSMU are linked to the foraging behaviour of the penguins for which performance data area available i.e. breeding, adult pygoscelids.
The authors cite Hinke et al. (2017) as the evidence supporting usage of the two gSSMU as appropriate strata. Pygoscelid penguins exhibit staggered breeding, with Adélies commencing first, followed by chinstraps then gentoos (Black, 2016). Adélie penguins are the first to fledge their chicks and thus cease to be centrally foraging, typically departing mid-February for their moulting grounds on the sea ice. Chinstrap penguins depart for a pre-moult foraging trip towards the end of February and return to land in order to moult, before departing again for their overwinter trip (Hinke et al. (2015); Hinke et al. (2019); Figure 2). Conversely, Gentoo penguins appear to remain in close proximity to their breeding colonies overwinter .
We use the Argos-CLS PTT telemetry data provided with this supporting study to characterise the actual at-sea habitat used, in the context of the relative stage of breeding for each species. For each species, we exclude locations with a “Z” error class and calculate the 99% Minimum Convex Polygon (home range) and associated area in \(nm^2\). For chinstrap penguins at Cape Shireff, this equates to a home range area of ~\(4,782nm^2\), or only 23% of the gSSMU to which their performance metrics are indexed against (Watters et al., 2020). For the same species at Copacobana the 99% MCP home range is 2,905\(nm^2\), or ~19% of gSSMU 1 in the Bransfield Strait. Similarly for Adélie penguins, the breeding foraging range occupied 1,139\(nm^2\) or only ~7% of the area of gSSMU #1. After breeding, available overwinter PTT telemetry and light geolocating data on chinstrap and adelie penguins suggests a wide dispersal westwards into the Pacific sector of the Southern Ocean, and eastwards into the Weddell Sea and Atlantic sectors, with a relatively small proportion of chinstraps from the study sites remaining within 500km of their breeding colonies (Hinke et al., 2019). Yet despite the evidence supporting widescale post-breeding migration of both Adélie and Chinstrap penguins, the model used by Watters et al. (2020) constrains both species from Copacobana to gSSMU 1 and Chinstraps from Cape Shireff to gSSMU 2 over winter (Watters et al. (2020); Supplementary Material 1 & 2, code lines 258 to 259). This has the effect of constraining the variability in performance indices from these species to LHR, LKB and ONI over winter in areas where the species has a demonstrated tendency to migrate out of. This is particularly important given that the fishery can now be characterised with a late autumn/early winter start (Figure 4).
Our preliminary review thus raises two areas of concern. Firstly, that the scales at which “local” predictors are summarised are in some cases five times greater than the habitat exploited by penguins. Secondly, that the known overwinter migratory behaviour of Adélie and Chinstrap penguins are poorly reflected in the model formulation. To demonstrate the impact that these ecological assumptions have on the model output, we rerun the model of Watters et al. (2020) with modified code. To avoid an overly burdensome paper, we shortly summarise those code changes here, and if requested during the meeting we are happy to include the rmarkdown version of this paper with the modified code in place.
However, we also note an additional coding error that may influence how the original (i.e. unmodified) results are interpreted. In summarising the model outputs into boxplots, the code relating to developing Figure 2 (Supplementary Material 1, lines 661-665) seemingly classifies the “Worst Case” (\({-0.5}\) \(^{\circ}\)C < ONI < 0.5 \(^{\circ}\)C; LKB > 1 Mt; and LHR \(\geqslant\) 0.1) using column 36 from the output dataframe, which actually reflects “ONI > 0.5 \(^{\circ}\)C; LKB > 1 Mt; and LHR \(\geqslant\) 0.1”; that is, ONI is in its “warm” phase as opposed to “neutral”. We thus relabel the output boxplot axis labelling to reflect the conditions representing “worst case” accordingly.
modifications
NEEDS WORKING ON - A LOT :) BULLET POINTS FROM THE GOOGLE DOC BELOW:
Our preliminary review of the evidence supporting localised effects of fishing coupled with broad-scale climactic phenomena having an impact on the vital statistics of pygoscelid penguins (performance and demographic trends) are based on assumptions that potentially do not reflect current knowledge of penguin breeding phenology and movement.
DESCRIPTION OF x6 SCENARIOS RUN WITH STIG
Of greatest concern, however, is that the interpretation of model outputs from both approaches (either from the original studies or the modified parameters we describe) are under boundary conditions that we feel are not appropriate. Both approaches consider only the fishery and broad-scale climate phenomena as the only two causes of krill abundance variability at geographic scales relevant to penguins. Neither study considers, for example, the impact of rebounding baleen whale populations or migratory male Antarctic fur seals beyond brief mentioning. Both taxa have increased in abundance throughout the life of the krill fishery, and there are sufficient telemetry and distance sampling studies in the scientific literature to demonstrate the degree and significance of spatiotemporal overlap with breeding penguin populations (see Santora and Veit (2013), Lowther et al. (2020), and Oosthuizen et al., Johannessen et al. and Lowther et al. submitted to this meeting, and Figure 3 as examples). Importantly, the distribution of these and numerous other unconsidered competitors is not uniform in either space or time and their impact on local availability of krill is likely to be considerable.
Similarly, the utilisation of broad-scale climatological phenomena to characterise impacts at scales that predators are dependent upon is problematic. The Amundsen Sea Low (ASL) is the dominant climate feature for the western Antarctic Peninsula. The El Niño Southern Oscillation (ENSO) modulates the the ASL, with El Niño (La Niña) shallowing (deepening) its pressure, causing more northwesterly (southeasterly) winds and upwelling (restricted influx) of Circumpolar Deep Water onto the shelf. The Southern Annular Mode also influences the pressure of the ASL, with the current trend of negative SAM constructively (destructively) interfering with ASL when in phase with El Niño (La Niña) events (e.g. Clem et al. (2016)). The result is a set of above-surface climate conditions that drive changes in water mass intrusion that are dependent on interactions between two climate processes. However the bathymetry of the Antarctic Peninsula is complex (particularly at scales that are important to centrally-foraging predators such as penguins) and the structuring of krill aggregations in time and space in the WAP have been linked to mesoscale circulation processes (Santora et al., 2012), which are unlikely to be uniformly affected by macroscale processes.
Our work into the future will progress along three lines, and we welcome any and all offers of collaboration into this work. Firstly, we will progress this debate into the scientific literature in order to ensure a balanced discussion occurs in that forum. Secondly, we will be examining in further detail some of the additional predictors used and their efficacy, the modelling frameworks into which they are brought, and how their incorporation influences the interpretation of the responses. Finally, we shall also be exploring alternative modelling approaches that reflect more of the physical and biological complexity of the system in question. In all cases, our goal is to ensure that the best available objective scientific evidence is presented to our environmental managers and, where appropriate, flag that disagreement exists. Our paper should be viewed in this light to generate constructive dialogue that addresses our common concern of the potential for localised fishing to impact dependent aspects of the ecosystem.
Penguin foraging behaviour derived from available ARGOS-CLS PTT data. A) Chinstrap penguins truncated at \(10^{th}\) March in line with known phenology (elongated grey track represents a single animal) B) Adélie penguins truncated to the end of January and C) Gentoo penguins until ~August, in line with known overwinter foraging behaviour. The SSMU are combined and coloured according to gSSMU (red; gSSMU 2, purple; gSSMU 1) with Chinstrap and Adélie penguin 99% MCP home ranges occupying between 7-19% of the gSSMU to which they were assigned.